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In monolayer hexagonal lattices, two inequivalent valleys appear in the Brillouin zone. With 
inversion symmetry breaking, we find chiral phonons with valley contrasting circular polarization 
and ionic magnetic moment. At valley centers, there is a three-fold rotational symmetry endowing 
phonons with a quantized pseudo angular momentum, which includes spin and orbital parts. From 
conservation of the pseudo angular momentum, crystal momentum and energy, selection rules in 
intervalley scattering of electrons by phonons are obtained. The chiral valley phonons are verified 
and the selection rules are predicted in monolayer Molybdenum disulfide. Due to valley contrasting 
phonon Berry curvature, one can also detect a valley phonon Hall effect. The valley-contrasting 
chiral phonon, together with phonon circular polarization, ionic magnetic moment, phonon pseudo 
angular momentum, valley phonon Hall effect, will form the basis for valley-based electronics and 
phononics applications in the future. 


Due to inversion symmetry breaking, valley¬ 
contrasting electronic physics has been proposed 


by Xiao et al. in 2007 [1 
dependent optoelectronics 


and followed by valley- 
and coupled spin and 


valley physics in monolayers of group-VI dichalcogenides 
AB 2 (A= Mo, W; B=S, Se)) lattices @|- Besides charge 
and spin, separated valleys in momentum space consti¬ 
tute another discrete degrees of freedom for electrons 
with long relaxation time, which leads to emergence 
of valleytronics, such as valley polarization and valley 
coherence on transition-metal dichalcogenides monolayer 
AB 2 materials 0 F 3 . The valley electron interband 
scattering involves a polarized photoexcitation and 
photoluminescence; however, the intervalley electron 
scattering will involve If-valley phonons [?]. Given the 
fact that electrons have definite chirality at valleys, a 
natural question then arises: whether do valley phonons 
have chirality and how does the chirality play a role in 
electronic intervalley scattering? 

In this Letter, we observe phonon chirality at two val¬ 
leys in honeycomb systems with broken inversion symme¬ 
try. By introducing isotopic doping [ll| or staggered sub¬ 
lattice potential [121 ] in graphene, the inversion symmetry 
can be broken. The monolayer hexagonal Boron Nitride 
which can be successfully grown over a wide area 0 , can 
also be expected to present valley chiral phonons. Besides 
these kinds of monolayer honeycomb AB lattices, mono- 
layer hexagonal AB 2 lattices which have been widely in¬ 
vestigated in valleytronics can also serve as paradigms 
for applications of K-v alley chiral phonons. The chi¬ 
rality of phonons at If-valleys not only decides the se¬ 
lection rules in electronic intervalley scattering but also 
can endow phononics with other potential effects, e.g., 
valley phonon Berry curvature and valley phonon Hall 
effect. And near if-valley centers phonon has extremes 
in dispersion and thus a large density of states. Therefore 
chiral valley phonons will play an important role in val¬ 
leytronics, especially in intervalley scattering of electrons 
or holes. 


Chirality of If-valley phonons. To study K valley 
phonons, we first focus on a two-dimensional honeycomb 
lattice model, where each unit cell has two sublattices A 
and B. The honeycomb AB lattice can serve as a sim¬ 
plified model to demonstrate general features of chiral 
phonons and can be generalized to other honeycomb ma¬ 
terials with broken inversion symmetry. The dynamical 
matrix is obtained by considering nearest neighbor inter¬ 
action. The spatial inversion symmetry is broken if the 
two sublattices A and B have different mass when a gap 
is open between the acoustic bands and optical bands as 
shown in the if and if' points in Fig. [Ha). From cal¬ 
culated eigenvectors, we can plot sublattice vibrations at 
valleys as shown in the insets of Fig. Oja). At valleys, all 
the vibrations are circularly polarized. For valley phonon 
modes of band 1 and 4 , both sublattices do opposite cir¬ 
cular motion; for band 2 and 3 while one sublattice is still 
the other sublattice does a circular motion. From valley 
if to if', all the circular motions will change to opposite 
directions. In latter discussion, we will focus on the if 
point while results at if' are achieved by time reversal 
symmetry. 

The phonon chirality can be characterized by polar¬ 
ization of phonon, which comes from the circular po¬ 
larization of sublattices. To consider the polarization 
along z direction, we look at the phonon eigenvectors 
e = (xi yi X2 y2) T (here we take a two-sublattice unit cell 
as an example, for a general case see Supplementary in¬ 
formation Sec. I). By defining a new basis where one 
sublattice has right-handed or left-handed circular po¬ 
larization as |Ai) = 77= (1 * 0 0) T ; |Li) = (1 — i 0 0) T ; 

| R2) = 77= (0 0 1 i) T ; IL2) = 7)5(001 — i) T , the phonon 
eigenvector e can be represented as 

n 

e = e Rc l-^a) + e L a \L a ), (1) 

a—1 

where e Ra = (i? a | e) = -j= (x a - iy a ), eL a = (L a \ e) = 
77 = (x a + iy a )- Then the operator for phonon circular 
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FIG. 1: Valley phonons in a honeycomb AB lattice, (a) 

Phonon dispersion relation of a honeycomb AB lattice. The 
insets show phonon vibrations for sublattices A and B in one 
unit cell at K' ( k x = — |^, k y = 0) and K ( k x = | ^,k y = 0), 
numbers 1 to 4 denote four bands. The radii of circles de¬ 
note vibration amplitudes, phase and rotation direction are 
included, (b) Phase correlation of phonon non-local part for 
sublattice A (upper two panels) and sublattice B (lower two 
panels) at ( K') (left) and K (right), (c) Phonon pseudo an¬ 
gular momentum for bands 1 to 4 at valleys K' and K. Here, 
the longitudinal spring constant Kl = 1, the transverse one 
Kt = 0.25, and m A = 1, m B = 1-2. The primitive vectors 
are (a,0) and (a/2, \/3a/2). 

polarization along z direction can be defined as 

n 

S Z = J2 d^> ( R «\ ~ l L “) (2) 

CK=1 

the phonon circular polarization equals to 

n 

s; h = e'S°eh = J2(\tR<*\ 2 -\tL a \ 2 )h, (3) 

a=l 

here n = 2 for two-sublattice unit cells. The phonon 
circular polarization can have a value between zth since 
J2 I e R a 1 2 + I £ L a \ 2 = 1- The Sp h has the same form 

a 

with that of phonon angular momentum j k a along z 
direction (see Supplementary information Sec. I). The 
phonon circular polarization comes from those of sub¬ 
lattices in a unit cell. Each sublattice can have a 
circular polarization s z a which equals to e'S^eh with 
= \ R c) {R a | - | L a ) (L a \. In Fig. [2(a), at valley K , 
s A = 0 , sfj = —h for band 2 while s z A = Ti, s z B = 0 for 
band 3, where the phonon circular polarization happens 


to be quantized; for band 1 and 4, sublattice A (B) is 
left (right) circularly polarized with |s^| < 1 (|s^| < 1 ). 
Therefore, phonon circular polarization s ph at valleys 
can be nonzero. If the sublattices have nonzero effective 
charges e a , one also can expect nonzero ionic magnetic 

n 

moment M z = ^ s z a e a /(2m a ) at valleys. By intro- 

OL—l 

ducing a staggered sublattice onsite potential, the val¬ 
ley phonon angular momentum can also be observed (see 
Supplementary information Sec. II). 

Phonon pseudo angular momentum at valleys. In 

a honeycomb lattice, at high symmetry points T, K , K' 
phonons are invariant under a 3-fold discrete rotation 
about the direction ( z ) perpendicular to the lattice 
plane. Under the rotation, one can obtain #l(^-,z)uk = 
e ~ l2 f l P h Uk ' where is defined as the pseudo angular 
momentum of phonon with wave function Uk • The phase 
correlation of the wave function comes from two parts, 
one is from local (intra-cell) part ek.a, another is from the 
non-local (inter-cell) part e lR, k . Thus under a 3-fold ro¬ 
tation, at valleys one can obtain spin pseudo angular mo¬ 
mentum I s for the local part and orbital pseudo angular 
momentum 1° for the non-local part. The math of pseudo 
angular momentum under 3-fold rotator is 1 + 1 = — 1 
since the eigenvalues of rotation are ± 1 , 0 . 

Under a 3-fold rotation, both circular polarization \R a ) 
and | L a ) are eigenstates of the operator with 

pseudo angular momentum l s R = 1 and l s L = — 1 re¬ 
spectively. Therefore, in Fig. [2(a), at valley K, l A = 
— 1, l s B = 1 for band 1 and 4, l s B = —1 for band 2 and 
l A = 1 for band 3. Since a 3-fold rotation center can 
be at any sublattice in a unit cell, the vibration of the 
rotation-center sublattice must be an eigenstate of the ro¬ 
tation operator 5R(^, z) due to no phase change for the 
orbital part of it. Therefore for non-degenerate modes, 
sublattices must do circularly polarized vibration other¬ 
wise they are still, which is consistent to our findings in 

Fig. [2 (a). 

The orbital pseudo angular momentum can be ob¬ 
tained from phase change under a 3-fold rotation, which 
is shown in Fig. [2 (b). We can obtain I a = t and 
l B = — r, r = ±1 labels the two valleys K and K' . 
Since phonon state must be an eigenstate of the rotation 
operator, pseudo angular momentum of phonon equals 
Iph = l A + l A = 1% + T both sublattices are vibrating; 
if one is still, it is decided by the other vibrating sub¬ 
lattice. Therefore we can obtain phonon pseudo angular 
momentum as listed in Fig. [2 (c). 

The chiral valley phonons can be observed in graphene 
systems with introducing isotope doping (please see Sup¬ 
plementary information Sec. Ill for detailed calculation), 
where the three possibilities of phonon spins are verified. 
The out-of-plane phonon modes will not play a role in 
phonon polarization along the direction perpendicular to 
the lattice plane. 
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FIG. 2: FF-valley phonons emitted in hole intervalley 
scattering in M 0 S 2 . The exciton is excited by a right polar¬ 
ized photon with energy (A) at K valley, (a) By absorbing a 
stimulated right-handed photon with energy (A kv + ha;), the 
excited hole in valence band is scattered to the other valley K' 
by emitting a FF-valley phonon with energy Tllo and pseudo 
angular momentum l p h = — 1. (b) A stimulated left-handed 
photon is absorbed and a phonon with l p h = 1 is emitted. The 
pseudo angular momenta of electrons in conduction band and 
valence band (l v ( c ) = ±1, 0) are also marked. 


Selection rules. It is well known that in graphene the 
double resonance D-peak in Raman spectrum is related 
to phonon modes in the vicinity of the K point [3, Eli . 
For electrons with zero moment along the normal direc¬ 
tion of the plane (z), the state is invariant under a 3- 
fold discrete rotation, the pseudo angular momentum is 
decided by the orbits on sublattice A or B. Based on 
the lattice structure in Fig. |T] (b), if we assume the va¬ 
lence band corresponds to the orbit on sublattice A and 
the conduction band corresponds to the orbit on sublat¬ 
tice B, we can obtain all the pseudo angular momenta 
^c(jj) = -F r - Therefore one can expect an azimuthal selec¬ 
tion rule l c — l v = At for interband transition by photons 
with right (<t+) or left (it—) circular polarization. 

A photon excited electron in one valley can decay to 
another valley by emitting a valley phonon through elec¬ 
tron phonon interaction. Due to circularly polarized mo¬ 
tion of sublattices at valleys, one can obtain a nonzero 
ionic magnetic moment, which can provide a possible 
mechanism for the electron phonon interaction. In the 
intervalley scattering by phonon the whole system has 
3-fold discrete rotation symmetry, thus we can expect 
a selection rule from the conservation of pseudo angular 
momentum, that is l c t v ){K)—l c i v \{K’) = ±1 by emitting 
a circularly polarized valley phonon (Z p h = ± 1 ), where 
momentum and energy conservations are also applied. 
Since 1 + 1 = -1 and K + K = K' (or K' + K' = K), a 
phonon with a pseudo angular momentum equal to that 
of the electron in the same valley. Thus we can expect 
an valley phonon with a specific pseudo angular momen¬ 
tum can be created. Due to the pseudo angular momen¬ 
tum of valley phonons are different, we can observe a 
circularly polarized infrared spectrum during the valley 


phonon interband scattering (please see Sec. IV of the 
Supplementary information for the detailed discussion). 

Due to spin-orbit coupling, the band structure of 
transition-metal dichalcogenides will have spin splitting 
A k at K valleys, e.g. for highest valence band at val¬ 
leys M 0 S 2 has a 150 meV splitting, while for the low¬ 
est conductance band only a negligible 3 meV splitting 
[3 . Thus we need consider the splitting A k for en¬ 
ergy conserving in electron intervalley scattering where 
the electronic spin keeps fixed. Furthermore, with a pho¬ 
ton absorbtion or emission we can expect to observe in¬ 
tervalley electron scattering at valley centers involving 
a valley center phonon, where the selection rules imply 
AZgi = A^ph + ^photon and Xk = A^-^ph A ?za^photonj where 
’+’ means emission and ’ means absorbtion. 

Recently experiments on valley polarization have been 
reported on some monolayer transition-metal dichalco¬ 
genides by optical pumping. In monolayer M 0 S 2 , Cao, 
et al. found that the pseudo angular momenta of elec¬ 
trons for valence and conduction bands are quite differ¬ 
ent from those in gapped graphene. They found valleys 
in the conduction band have Al pseudo angular momen¬ 
tum, while both valleys in the valence band have 0 as 
shown in Fig. [2] thus they found the absorption of circu¬ 
larly polarized photons at valleys. 

In Fig.[2l with spin-orbit coupling M 0 S 2 has a bandgap 
of A = 1.65eV and a spin-splitting for the highest valence 
band at K valley A Kv = 150 meV 0 . Through a right- 
handed polarized photon a pair of exciton are excited 
at K valley, where the blue lines correspond to the ab¬ 
sorption of a right-handed photon with energy A. Since 
the excited electron is in the valley center, which can¬ 
not be scattered to another valley through emitting a 
phonon. However, due to the large spin-splitting of va¬ 
lence band the hole can be scattered to another valley 
by a circularly polarized photon and phonon by keeping 
spin fixed. Using the Quantum-Espresso code 0, we 
obtain phonons for all bands at valleys as shown in Ta¬ 
ble U (please see Supplementary information Sec. V for 
detailed calculation). At K valley, phonons with energy 
14.4 and 48.0 meV have pseudo angular momentum of 
Z p h = —1 while the mode 48.0 meV phonon has a lit¬ 
tle polarization, and phonons with energy 21.5, 40.0 and 
45.1 meV have pseudo angular momentum of l p h = 1. 

Through a right-handed light scanning on the sample, 
we can observe some resonance peaks on A Kv A hut p h 
where stimulated photons are emitted, as shown in Fig. 
[2])a) . The emitted phonon will have a pseudo angular 
momentum of Z p h = — 1 at K point due to the selection 
rules at valleys. Thus we can observe a resonance peak 
of 164.4 meV (another peak of 198.0 meV is not obvi¬ 
ous due to the small polarization of phonon). And for a 
stimulated left-handed photon in Fig. [2])b) , we can only 
observe one peak of 190.0 meV correspond to phonon 
mode of 40.0 meV while the other two modes will not be 
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TABLE I: Chiral phonons in K Valley of M 0 S 2 . Mono- 
layer M 0 S 2 have 9 modes (n) of phonon with energy huj 
(meV), Mo and S are located in A and B respectively as show 
in Fig. 1 (b) in the main text, thus the orbital pseudo angu¬ 
lar momentum Im 0 = 1 and lg = —1. Sm d (ss)j So (£s), and 
l p h are circular polarizations of Mo (S), spin pseudo angular 
momenta of Mo (S), and phonon pseudo angular momentum. 
The Mirror symmetry (Ms) is relative to the plane of the 
monolayer of M 0 S 2 ; it is 1 (-1) if the mode is even (odd) 
under the mirror symmetry operation. 


n 

hoj 

s Mo 

S S 

1 S 

l Mo 

I s 

L S 

^ph 

M s 

1 

14.4 

0.64 

0 

1 

0 

-1 

1 

2 

21.5 

0 

- 0.34 

0 

-1 

1 

-1 

3 

31.9 

-0.18 

0.41 

-1 

1 

0 

1 

4 

40.0 

0 

-0.50 

0 

-1 

1 

1 

5 

40.2 

0 

0 

- 

- 

- 

-1 

6 

41.5 

0 

0.50 

0 

1 

0 

1 

7 

45.1 

0 

-0.40 

0 

-1 

1 

-1 

8 

48.0 

0.06 

0 

1 

0 

-1 

1 

9 

49.2 

-0.34 

0.33 

-1 

1 

0 

1 



FIG. 3: Phonon Berry curvature and valley phonon 
Hall effect in a honeycomb lattice, (a) Berry curvature 
of band 1 (bottom contour plot) and band 2 (top 3D plot). 
Band 3 (Band 4) has phonon angular momentum opposite to 
that of band 2 (band 1). (b) ((c)) Schematic of valley phonon 
Flail effect (Hall current denoted by olive curve arrows) under 
a strain gradient (orange arrows), where valley phonons are 
excited by a ray of right-handed or left-handed polarized light 
(red wave lines). The parameters are the same with those in 
Fig. 1. 


involved since they are odd under mirror operation (see 
Table I) while in the electron-phonon scattering process 
the electron wave function keeps even under a mirror 
operation relative to the x-y plane. During these two 
processes, only iC-valley phonon can be emitted while a 
chiral phonon at K' can be emitted in the scattering of 
hole of an exciton at K' valley which can be excited by 
a left-handed photon with 1.65 eV. After scattered by 
phonon and photon, the pair of electron and hole locate 
in different valleys and can have a long life time; its con¬ 
densation has been theoretically observed very recently 
[20]. Therefore, a specific chiral phonon at a definite val¬ 
ley can be obtained through a stimulated photon. With 
the two-step polarized light shinning on the sample, a 
large number of valley phonons with definite frequencies 
can be created. 

Valley phonon Hall effect. In the presence of an 
in-plane electric field, an electron will acquire an anoma¬ 
lous velocity proportional to the Berry curvature in the 
transverse direction 2l|, |22|. Recently the electronic 
valley Hall effect proposed in jT|| has been experimen¬ 
tally observed in monolayer M 0 S 2 transistors [23( and in 
graphene superlattices [24|. As discussed above, phonons 
with definite frequencies at a specific valley can be mas¬ 
sively created, thus for valley phonons, if its Berry cur¬ 
vature is nonzero we can also expect to observe valley 
phonon Hall effect in the presence of an in-plane gradient 
strain field. Such valley phonon Hall effect can provide 
us another way to observe valley phonons. 

With the breaking of spatial inversion symmetry we 
observe nonzero phonon Berry curvature at valleys as 
shown in Fig. [3] (a) (see Supplementary information 
Sec. VI for derivation). Band 2 and band 3 have large 


Berry curvatures at valleys, while those of band 1 and 4 
are small. Due to the nonzero phonon Berry curvature, 
applying a strain gradient Strain along x direction, 
phonons excited at a different valley will go to a different 
transverse direction since u anom oc —^strain x f2 in anal¬ 
ogy to electrons. If the photon polarization is reversed, 
the transverse phonon current would be reversed as 
shown in Fig. [3] (b) and (c). 
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